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Abstract — A multi-hop relay channel with multiple antenna 
terminals in a quasi-static slow fading environment is considered. 
For both full-duplex and half-duplex relays the fundamental 
diversity-multiplexing tradeoff (DMT) is analyzed. It is shown 
that, while decode-and-forward (DF) relaying achieves the opti- 
mal DMT in the full-duplex relay scenario, the dynamic decode- 
and-forward (DDF) protocol is needed to achieve the optimal 
DMT if the relay is constrained to half-duplex operation. For 
the latter case, static protocols are considered as well, and the 
corresponding achievable DMT performance is characterized. 

I. Introduction 

Relays are commonly used in wireless networks to improve 
performance, although the fundamental capacity limits of relay 
channels have yet to be fully characterized, even for simple 
systems [1]. Rather than focus on capacity limits, we are 
interested in characterizing the tradeoff between rate gain 
through multiplexing versus the robustness gain through di- 
versity associated with multiple-antenna relays. We will focus 
on a multiple antenna multi-hop system in which the source 
transmission can only be received by the relay terminal, as 
shown in Fig.[T] We call this the multiple-input multiple-output 
(MIMO) multi-hop relay channel. The links are assumed to 
be quasi-static, frequency non-selective Rayleigh fading, and 
the channel state information (CSI) is available only at the 
receiving end of each transmission. 

We analyze this system in terms of the diversity- 
multiplexing tradeoff (DMT) in the high signal-to-noise ratio 
(SNR) regime introduced in [2]. DMT analysis is useful in 
characterizing the fundamental tradeoff between the reliability 
and the number of degrees of freedom of a communication 
system. In DMT analysis, reliability is measured in terms of 
the diversity gain, which characterizes the rate of decay of the 
error probability with increasing SNR. The degrees of freedom 
is measured by the spatial multiplexing gain, which is the 
rate of increase in the transmission rate with SNR. While the 
DMT analysis is a tool to characterize the fundamental limits 
of a communication system in a fading environment, practical 
space-time codes that approach these theoretical limits have 
been designed [3]- [6]. 

In a channel with relays, the source's transmission is re- 
ceived by both the relays and the destination, and the source 
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Fig. 1. The (Mi, M2, M3) MIMO multi-hop relay channel model consid- 
ered in the paper. There is no direct link from the source terminal (S) to the 
destination terminal (D). 



and the relay terminals cooperate to transmit the message 
to its intended destination [7], [8]. DMT analysis has been 
extensively applied to this general relay channel model; how- 
ever, a full characterization of the DMT curve is still an open 
problem. In [9] the DMT of the half-duplex single-antenna 
relay channel is analyzed and a dynamic decode-and-forward 
(DDF) protocol is proposed. In DDF, the relay terminal listens 
to the source transmission until it can decode the message, and 
then starts transmitting the message jointly with the source 
terminal. The DMT of DDF is shown to dominate that of all 
other protocols, but for high multiplexing gains it does not 
meet the cut-set upper bound, which dictates the maximum 
possible of such gains [10]. In [11] DDF performance is 
improved slightly by using superposition coding. In [12], under 
the assumption of full CSI at the relay terminal, the compress- 
and-forward protocol is shown to achieve the optimal DMT 
performance. There has also been some recent interest in 
the DMT analysis for multi-hop relay systems; in [13] and 
[14] multiple single antenna relays operating in a distributed 
manner are considered. Due to the distributed nature of the 
relay nodes, amplify-and-forward relaying is considered, under 
which the achievable DMT is characterized. 

In contrast to these prior works, we consider a MIMO multi- 
hop relay channel. For this model, the relay can decode the 
message without sacrificing degrees of freedom. While we 
only derive results for a single relay in this paper, our results 
can be extended to multiple relays. In the case of a full-duplex 
relay, we show that the DF protocol with a block Markov 
structure achieves the optimal DMT performance. In the half- 
duplex relay case, we first find the DMT of static protocols 
in which the source and the relay transmission periods are 
fixed, independent of the channel realization. On the other 



hand, we show that the DDF protocol of [9], in which the 
time allocation depends on the realization of the source-relay 
channel, achieves the optimal DMT performance. In the multi- 
hop scenario, since the relay and the source do not transmit 
simultaneously, they do not need to use distributed space- 
time codes, which are harder to realize in practice [15], [16]. 
Furthermore, there is no need to inform the source or the 
destination terminals about the relay decision time as opposed 
to the general relay scenario. Hence, the dynamic relaying 
scheme in the case of the multi-hop relay channel can be 
realized by using an incremental redundancy code at the 
source [6], and any DMT-optimal space-time code at the relay. 
Although the DMT of DDF has been previously shown to 
dominate that of other protocols in the case of general half- 
duplex relay channels, our results prove its optimality in the 
multi-hop relay scenario. In a concurrent work [17], Gharan 
et al. prove the optimality of the DDF protocol in a single- 
antenna multiple access relay network. 

II. System Model 

We consider a three node multi-hop channel composed of 
source, relay and destination terminals with Mi, M 2 and M 3 
antennas, respectively, as in Fig. Q] We call this system an 
(Mi, Mi-, M3) multi-hop relay channel. The source-relay and 
the relay-destination channels are given by 

Yj = y^HiXi+Wi, (1) 

for i = 1, 2, respectively, where Yj, i — 1,2, are the received 
signals at the relay and the destination, respectively. Note that 
the source transmission is not received at the destination in 
our multi-hop relay channel model. Channels are assumed to 
be frequency non-selective, quasi-static Rayleigh fading and 
independent of each other; that is, for i = 1,2, H; is an 
Mj+i x Mi channel matrix whose entries are independent 
and identically distributed (i.i.d.) complex Gaussian random 
variables with zero means and unit variances (i.e., they are 
CN(0, 1)). The additive white Gaussian terms also have i.i.d. 
entries with CAf (0,1). Xj, i = 1,2, are Mi x T source 
and relay input matrices, where T is the total number of 
transmissions over which the channel is constant. We have 
short-term power constraints at the source and the relay given 
by trace(E[Xf Xi]) < M{T. For % = 1,2, we define 

M* = min{Mj, M l+ i}. 

We assume that only the receivers have channel state infor- 
mation. 

Following [2], for increasing SNR we consider a family of 
codes and say that the system achieves a multiplexing gain of 
r if the rate R(SNR) satisfies 

R(SNR) 

lllTl = T 

SNR^oo \og(SNR) 

The diversity gain d of this family is defined as 

\o g P e (SNR) 
SNB^oo log(SNR) ' 



in which P e (SNR) is the error probability. For each r, define 
d(r) as the supremum of the diversity gain over all families 
of codes. The full characterization of the DMT curve for a 
MIMO system is given in the following theorem [2]. 

Theorem 2.1: For a MIMO system with M\ transmit and 
M2 receive antennas and sufficiently long codewords, the 
optimal DMT curve dM 1 ,M 2 ( r ) is given by the piecewise- 
linear function connecting the points (k,d(k)), k = 
0, . . .,min(Mi,M 2 ), where d(k) = (M x - k)(M 2 - k). 

For the rest of the paper, we always consider codes with 
sufficiently long codewords so that the error event is dominated 
by the outage event. 

III. DMT of MIMO Multi-hop Relay channels 

A. Full-duplex Relaying 

We first consider the full-duplex relay case. The next 
theorem shows that the DMT tradeoff of the end-to-end system 
is equal to the worst-case DMT tradeoff of each link along 
the multi-hop path. The DMT characterization given here for 
a single relay can be easily generalized to multiple full-duplex 
relays. 

Theorem 3.1: The DMT < ^L 1 M 2 M 3 ( r ) °^ an 
(Mi, M 2 , M3) full-duplex system is characterized by 

d M 1 ,M 2 ,M 3 ( r ) = niin{d A/l ,M 2 (r),dM 2 ,Af3(^)}- ( 2 ) 
Proof: The result follows easily as DF achieves the 
capacity for a full-duplex multi-hop relay channel [1]. ■ 

B. Static Protocols for Half-duplex Relaying 

In the half-duplex relay scenario, the total T time units need 
to be divided among the source and the relay transmissions. 
We first consider static protocols for which the time allocation 
is fixed, independent of the channel states. However, similar 
to the generalized decode-and-forward protocol in [18], we 
consider unequal division of the time slot among the source 
and the relay. The source transmits during the first aT channel 
uses, where < a < 1. The relay tries to decode the message 
and forwards over the remaining (1 — a)T channel uses. We 
call this protocol decode-and-forward with fixed allocation 
(fDF), and its DMT is given in the next proposition. 

Proposition 3.2: The DMT of the half-duplex 
(Mi,M2, M3) relay channel with fixed time allocation 
a (0 < a < 1) is 

d M?M 2 M 3 ( r ) = min jrf Ml ,M 2 (^) ,d M2 ,M 3 (j~T^) } (3) 
Proof: This result follows easily from Theorem 12. II with 
simple scaling of the DMT curve due to time division. ■ 

We can see from the above DMT that the highest multi- 
plexing gain for the fDF scheme is min{aM^, (1 — a)M^}. 
On the other hand, the highest diversity gain is limited to 
M 2 min{Mi, M 3 }. We illustrate the DMT of a (4, 2, 3) system 
with a fixed time allocation of a = 0.3 in Fig. [2] 

Since different time allocations result in different DMT 
curves, we can optimize the time allocation based on the 
multiplexing gain [18], [19]. We call this protocol DF with 
variable time allocation (vDF). Note that this is still a static 
protocol since the time allocation variable is determined based 
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Fig. 2. The dotted and the dashed curves correspond to d4,2(r) and efo,3(r), 
respectively. Note that the dashed curve also corresponds to the DMT in the 
case of a full-duplex relay terminal. The solid curve is the DMT curve of a 
(4, 2, 3) half-duplex multi-hop relay with the fDF protocol and a = 0.3. 



Fig. 3. The DMT curve of a (4, 1, 3) multi-hop relay channel. The two 
topmost curves correspond to the cut-set bounds, where the dashed curve is 
also the DMT for a full-duplex relay. The DDF, vDF and fDF protocol with 
a = 0.5 are also illustrated, where the DDF curve is the optimal DMT with 
half-duplex relaying. 



only on the multiplexing gain and is independent of the 
channel realization. For each multiplexing gain r, the diversity 
gain is the minimum of the two diversity gains in (O; hence the 
optimal time allocation variable a(r) is the one that satisfies 



jvDF 

a M u M 2 ,M 3 



(r) 



1 M 1 ,M 2 



a{r) J \ 1 — a[r) 

Corollary 3.3: The number of degrees of freedom of 
an (Mi,M 2 ,M 3 ) multi-hop relay channel with the vDF 
protocol is m*+M' 1 wr, il e the maximal diversity gain is 
M 2 min{Mi,M 3 }. 2 

We now present the DMT for some special cases because 
a general closed form expression is not tractable. We first 
consider the (Mi,l,Mg) system. Since the hops for this 
setup are multiple-input single-output (MISO) and single-input 
multiple-output (SIMO) systems, the DMTs are characterized 
as d,Mi,M i+ i — M£(l — r), i = 1,2. From (HJi and defining 
A = M x */M 3 * and B = 1 — r - A(l + r), we find 



(4) 



a{r) 



-B + yjB" 1 - AA(A - 1) 



2(A-l) 

for A j= 1. We have a(r) = 0.5 if A = 1. The DMT achieved 
by the vDF protocol in a (4, 1, 3) system is plotted in Fig. [3] 
In this figure, we also plot the DMT for the fDF scheme with 
a fixed time allocation a = 0.5. 

If we have Mi = M3 = M, then the optimal time allocation 
is a — 0.5 independent of the multiplexing gain, and the DMT 
is given by d v ^^ M (r) = d M ,M 2 (2r). 

C. Dynamic Decode-and-Forward Protocol for Half-duplex 
Relaying 

In [9], Azarian et al. proposed the dynamic decode-and- 
forward protocol for the cooperative relay channel with single 



antennas. In DDF for the relay channel, the source transmits 
during the entire timeslot using an incremental redundancy 
type codebook. This code design enables the relay to decode 
the message after receiving only a portion of the codeword; 
hence the relay decodes the message when the accumulated 
mutual information over the source-relay channel is sufficient 
for the transmission rate. Thus, the relay decoding time 
becomes a random variable that depends on the source-relay 
channel quality. As soon as the relay decodes the message, it 
starts transmitting. 

The achievable DMT of the DDF scheme in the case of 
a single-antenna cooperative relay channel is characterized in 
[9], where it is shown to dominate the DMTs of amplify-and- 
forward (AF) and decode-and-forward (DF) based protocols 
and, more strikingly, to achieve the DMT upper bound for 
multiplexing gains r < 0.5. Hence, DDF is DMT-optimal in 
this range of low multiplexing gains for the single antenna 
cooperative relay channel. 

Here, we consider using the DDF protocol for the multi- 
antenna multi-hop relay channel, and show that it achieves 
the DMT cut-set upper bound; that is, DDF is DMT-optimal. 
The intuitive explanation behind the optimality of DDF in 
this setting is as follows: In the multi-hop relay scenario, 
the message needs to be decoded at the relay terminal, since 
otherwise the destination would not be able to decode it either, 
due to the data processing inequality. However, any fixed time 
allocation scheme either wastes multiplexing gain, since it 
cannot utilize the good states of the source-relay channel, or 
results in outage in the case of a poor quality (low SNR) 
source-relay channel. DDF, by enforcing decoding at the relay 
and dynamically allocating the source transmission time based 
on the source-relay channel state, achieves the optimal DMT 



performance. 

Theorem 3.4: For the (Ml, M 2 , M 3 ) system with rate R = 
r log SNR, the outage probability of the DDF is given by 

P out {r) = SNR~ dDDF ^) 

where 

2 M* 

d DZ?i » = inf V V(2j - 1 + |M - M i+1 |)a i;i 

(5) 

and 

2 4 j(a a ,a 2 )eft M *+xft M 2* + 
in which we have defined 

M* 

Siicxi) 4 5^(1 - Oij) + . for * = 1,2. (6) 

j=i 

Proof: The proof can be found in Appendix U ■ 

The optimality of the above DMT achieved by the DDF 
protocol is shown in the following theorem. 

Theorem 3.5: DDF is DMT-optimal for MIMO multi-hop 
half-duplex relay channels. 

Proof: The proof can be found in Appendix HI] ■ 

Corollary 3.6: The number of degrees of freedom of an 
(Mi,M2,M3) multi-hop relay channel is m'+m* ' wn il e the 
maximal diversity gain is M2 min{Mi, M3}. Hence, the end- 
points of the DMT curve can also be achieved by static 
relaying, i.e., with fixed time allocation corresponding to the 
multiplexing gain. 

It can be seen from Theorem 13.41 that the DMT of a 
half-duplex multi-hop relay channel is not a piecewise-linear 
function as in the case of a point-to-point MIMO channel. 
While it is hard to give a general closed form expression for 
the DMT of MIMO multi-hop channels, for given Mi,M 2 
and M3 and a fixed multiplexing gain r, the optimization 
problem in © can be converted into a convex optimization 
problem, and hence can be solved efficiently [20]. We now 
give an explicit characterization of the DMT for some classes 
of multi-hop relay channels. 

Corollary 3.7: The DMT of an (Mi, 1, M 3 ) system is 

dSf ) f, M3 (r)=min(Mi,M 3 )^^ 

1 — r 

for < r < 1/2, and elsewhere. 

In Fig. [3] we illustrate the DMT of the (4, 1,3) multi-hop 
relay channel, which is achieved by the DDF protocol. We 
see that the DDF DMT dominates that of the static protocols 
at all multiplexing gains except the end-points. As stated in 
Corollary 13.61 these end-points can be achieved by the static 
fixed time allocation fDF protocol as well. 



DMT of (2,2,2) half-duplex relay channel 
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Fig. 4. The DMT of a (2, 2, 2) system. From top to bottom, the three 
curves correspond to the full-duplex relay DMT, the half-duplex relay DMT 
which is achievable by DDF protocol, and the DMT of the static protocol 
with a = 0.5. 

Corollary 3.8: The DMT of the (2, 2, 2) system with a half- 
duplex relay is given by 

f if < < 1/2 

d DDF {r) = i if 1/2 < b < 2/3 (7) 

[ 4^_r) if 2 /3<6<l 
The DMT of the (2, 2*2) system is plotted in Fig. g] The 
topmost curve in the figure is the DMT of a 2 x 2 MIMO 
system, which can be achieved by a full-duplex relay. The 
lowest curve is the DMT of the vDF protocol. Note that for 
this symmetric scenario vDF reduces to fDF with a = 0.5. 

IV. Conclusions 

We have derived the diversity-multiplexing tradeoff of 
MIMO multi-hop relay channels for both full-duplex and half- 
duplex relays. For full-duplex relays, it is easy to show that 
the decode-and-forward protocol achieves the optimal DMT, 
which is simply the minimum of the DMTs of the two links. 
This applies to multiple relays as well; that is, the DMT of 
the end-to-end system will be limited by link with the smallest 
DMT. In the case of a half-duplex relay, we have shown that 
the dynamic decode-and-forward protocol, in which the relay 
listens until decoding and then forwards, achieves the optimal 
DMT, which is no longer a piecewise-linear function of the 
multiplexing gain. We have also shown that this optimal DMT 
performance cannot be achieved by static time allocation. 
Finally, we have provided explicit expressions for the DMT 
of some classes of half-duplex multi-hop relay systems, and 
compared the achievable performances with fixed and dynamic 
time allocation. 

Appendix I 
Proof of Theorem I3.4I 

For the achievability scheme, we assume that the inputs at 
both the source and the relay are Gaussian with identity covari- 



ance matrices. Let the transmission rate be R = 7-logSNR, 
and define 

C< (H t ) 4 logdet (i + ™HiHt) • (8) 

The relay listens for aT channel uses until it decodes the 
message. Hence, we have 

rlogSNR 



If a > 1 then the relay is in outage, which leads to an outage 
for the whole system. If a < 1, then the relay transmits during 
the rest of the timeslot for (1 — a)T channel uses. Conditioned 
on successful decoding at the relay, i.e., a < 1, the outage 
probability over the second hop is given by 



P{rlogSNR> (l-o)C 2 (H 2 )} 

P<rlogSNR>[l-?|^)<7 2 (H 2 ) 



P{ rlogSNR > 



Ci(Hi)0 2 (H 2 ) 
C 1 (H 1 ) + C 2 (H 2 ) 



(9) 



Let Aj i, . . . , \i.m» be the nonzero eigenvalues of ILHj 
for i = 1, 2. Suppose X tJ = SNR~ a * J for j = 1, . . . , M*, 
i = 1,2. We haveQ 



C,(H,) = logJJll+^A^ 

3=1 V 

M* 

= log]JSNR (1 ^ )+ 

3=1 



(10) 



where (x) + = max{0,a;}. Using these exponential equalities, 
we can rewrite © as follows 



P{r log SNR> (l-o)C 2 (H 2 )} 

Cri(H!)C7a(H 2 ) 



P I rlogSNR > 
pj log SNR r > 
P{r> 



Ci(Hi) + C 2 (H a ). 
logSNR Sl(Ql) logSNR S2(Q2) 
logSNR Sl(ai) +logSNR S2(Q2) 
Si(a 1 )S 2 (a 2 ) 



S 1 (a 1 ) + S 2 {a 2 ) 



where we have Si(oti) = J2j=i(l — a hj) + - 

Then the overall outage probability can be written as 



Pout(r) = P{r>S' 1 (a 1 )} 

+p{Si(ai) >r> 



Si(a 1 )S 2 (cx 2 ) 
+5 2 (a 2 ) 



'Define the exponential equality as f(SNR)=SNR c , if 
limsjvfl—,00 '"log jfjvj^ = c - The exponential inequalities < and > 
are denned similarly. 



We define 

01 = {(ai,aa):r>5i(ai)} 



02 = |(ai,a 2 ) : Si(qi) > r > 

= 01 U 2 



5 , i(qi)5 2 (Q2) 
S , i(a 1 ) + S' 2 (Q !2 ) 



Then using the joint probability of the eigenvalues of HjH j 
given in [2], the outage probability can be computed as 

P out {r) = / p(a 1 ,a 2 )da 1 da 2 
Jo 

r 2 Mf 

= / nn sM " (2i " i+|Mi " Mi+iilay(iai(ia2 
»=ii=i 



where 0' = n (ft M *+,ft M =+). 

Using Laplace's method as in [2], we obtain the expo- 
nential behavior of the outage probability as P ou t(r) = 
SNR- dDDF( - r \ where 



d UUb {r)= inf /(ai,a a ) 

(cti,a 2 )eO' 



(11) 



and 



2 A/* 

/(ai,o a ) = Y^C*3 - 1 + |M< - M i+1 \)a id . (12) 

»=1 3=1 

Next, we define 

01 = {K,a 2 )eK M * + xK M Z + \ 

oh,i > ■ > OL it M? > 0,r > 5i(a x )} 

02 = {(ai,a a ) eft M ?+ x^ M2 * + | 



5i(ai) + 5 2 (a a ) 



a»,i > ■ ■ • > Q!i,Mf > 0, r > 
= 0i U 2 

We can see that 0' = 0. Hence, 

d DDF {r)= inf ./( ai ,a 2 ) 

(ai,o 2 )£(3 

= inf< inf /(a 1 ,a 3 ), inf /(ai,a 2 ) 

[(cn,a 2 )ee>i (ai,oa)EO a 

= inf /(ai,a 2 ) 

(ai,a 2 )ed 2 

in which the last equality follows since we have 

i>i(ai) + ^2(a2) 
for all (ai,a 2 ), and hence 0i C 2 . 



Appendix II 
Proof of Theorem |3.5I 

We first give an upper bound for the DMT of the MIMO 
multi-hop half duplex relay channel, and show that the DDF 
DMT given in Theorem 13.41 matches this upper bound. Let 
a G (0, 1] be the portion of the source transmit time, i.e., 
the source transmits over the first aT channel uses. Hence, 
the relay transmits over the remaining (1 — a)T channel uses. 
Here we assume that the time allocation is independent of the 
message, i.e., it cannot be used for information transmission. 
As shown in [12] this does not affect the DMT of the system. 

From the two cut-set bounds, the instantaneous capacity 
C(Hi,H2) is upper bounded by 

max mm{a7(X 1 ;Y 1 |H 1 ),(l-a)/(X 2 ;Y 2 |H 2 )}. 

a,Px. 1 Px 2 

Since the capacity is maximized with Gaussian inputs, and 
linear scaling of the power constraint does not affect the high 
SNR analysis, the instantaneous capacity can be bounded as 

C(Hx,H 2 ) < maxmm{a(5i(Hi), (1 - a)C 2 (H 2 )}, 

a 

where we define 

Q(Hi) 4 logdet(I + SNRHjHt). 

We can further upper bound the capacity by assuming optimal 
time allocation at each channel realization. The instantaneous 
capacity is maximized at each channel realization for 

«(H 1 ,H 2 )= - - ^ (H2 - } . 
V 21 Ci(Hi)+C 2 (H 2 )' 

and the corresponding upper bound is 

Ci(Hi)<5 2 (H 2 ) 



C(H X ,H 2 ) < 



Ci(Hi)+C 2 (H 2 ) 

For a transmission rate of R = r log SNR, the outage 
probability lower bound is given by 

C 1 (H 1 )C 2 (H 2 ) \ 
C 1 (H 1 )+C 2 (H 2 )/' 1 

Using the characterization of the eigenvalues of the channel 
matrices given in Appendix Q] we obtain 



Pout(r) > P \ r log SNR > 



Poutir) > P\r> 



Si(ai)S 2 (a 2 ) 
5i(ai) + S 2 (a 2 ) 
in which S,-(aj) is as defined before. Next, we define 

S 1 (a 1 )S 2 (cx 2 ) 



O 



(a, (3) : r > 



S 1 {a 1 ) + S 2 (a 2 
Then the outage probability is lower bounded by 



Pout(r) > 



p(a, /3)dad(3 



2 M* 



Similar to Appendix U we obtain the exponential behavior 
of the above outage probability using Laplace's method. Note 
that, since O' — C> 2 , the outage probability upper bound has 
the same diversity gain function as the DDF protocol. Hence, 
DDF is DMT-optimal. 
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